Molecular engineering of antiferromagnetic rings for quantum computation 
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The substitution of one metal ion in a Cr-based molecular ring with dominant antiferromagnetic 
couplings allows to engineer its level structure and ground-state degeneracy. Here we characterize a 
CrrNi molecular ring by means of low-temperature specific-heat and torque-magnetometry measure- 
ments, thus determining the microscopic parameters of the corresponding spin Hamiltonian. The 
energy spectrum and the suppression of the leakage-inducing S-mixing render the CryNi molecule 
a suitable candidate for the qubit implementation, as further substantiated by our quantum-gate 
simulations. 

PACS numbers: 75.50.Xx, 03.67.-a, 75.40.Cx 



Due to their relative decoupling from the environment 
and to the resulting robustness, electron spins in solid- 
state systems are currently considered among the most 
promising candidates for the storing and processing of 
quantum information (QIP) 0. In this perspective, 
an increasing interest has recently been attracted by a 
novel class of molecular magnets, including both ferro- 
magnetic and antiferromagnetic Q systems. In the 
latter case the quantum hardware is thought as a collec- 
tion (e.g., a planar array) of coupled molecules, each cor- 
responding to a different qubit. A major advantage with 
respect to analogous schemes based on single-spin encod- 
ings would arise from the larger dimensions of the physi- 
cal subsystem, and on the resulting reduction of the spa- 
tial resolution which is required for a selective addressing 
of each qubit by means of local magnetic fields fH. A 
detailed, system-specific investigation is however manda- 
tory in order to verify the actual feasibility of this ap- 
proach, and to suitably engineer the intra- and inter- 
cluster interactions, the coupling between the computa- 
tional and the environmental degrees of freedom, as well 
as the gating strategy. 

It is the purpose of the present Letter to argue the 
suitability of Cr-based antiferromagnetic molecular rings 
for the qubit implementation on the basis of a detailed 
theoretical and experimental investigation of its wave- 
functions and energy levels, and of its simulated time 
evolution as induced bysequences of pulsed magnetic 
fields. Molecular rings [J| are characterized by a cyclic 
shape and by a dominant antiferromagnetic coupling be- 
tween nearest neighbouring ions. In the absence of ap- 
plied fields and for even numbers of spin centers, their 
energy spectrum typically consists in a singlet ground 
state and in characteristic rotational excitations 5] . The 
recently demonstrated substitution of a Cr^+ ion with a 
divalent transition metal provides an extra spin to 



the otherwise fully compensated molecule: this may in 
turn result in the formation of a ground state doublet 
energetically separated from the higher energy levels, i.e. 
in a suitable level structure for the qubit implementation. 

In the octanuclear heterometallic ring of our present 
concern, one of the Cr'^+ (s=3/2) ions is substituted by 
a Ni^+ (s=l) one. The cyclic molecule, with a diameter 
d I nm, is characterized by a planar arrangement; we 
thus define 6 to be the angle between the static magnetic 
field Bo and z, the latter being perpendicular to the ring 
plane. The spin Hamiltonian corresponding to the single 
molecular magnet reads 

8 8 

H = ^ Ji Si • s,+i + ^ d» [si, - Si(s, -I- 1)/3] + 

i=l i=l 

8 8 

^ Si ■ Dy • Sj + MS X! ' 

i<j—l i—1 

where the first term accounts for the dominant isotropic 
exchange interaction, the second and third ones describe 
the anisotropic local crystal-field and the intracluster 
dipole-dipole interaction, respectively; isotropic g factors 
are assumed for the last, Zeeman term. The anisotropic 
part of 7i does not commute with the squared total spin 
operator S^, and thus mixes subspaces corresponding to 
different values of the total spin (S-mixing) . Due to its 
reduced symmetry, Ti. can no longer be independently di- 
agonalized within each (25" -t- l)-dimensional block: an 
efficient solution scheme, based on an irreducible tensor 
operator formalism, has therefore been developed (see 
Ref. [j] and references therein) and applied to the present 
case. 

In order to estimate the parameters entering the above 
spin Hamiltonian, we measure the heat capacity C as 
a function of the temperature (0.4 < T < 10 K) and 
of the magnetic field (0 < Bo < 7 T) 0. The sam- 
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FIG. 1: Specific heat C / R of a single CrTNi crystal as a func- 
tion of the temperature (T) and of the magnetic field (-Bo). 
We attribute the slight discrepancy between the theoretical 
curve and the experimental data at Bo = and low T to 
an experimental artifact, probably due to the presence of de- 
fected rings within the sample. Inset: torque signal measured 
on aO. 7x1x0.5 mm^ CryNi crystal at T = 2 K and Bo = 2 T. 
In our experimental set-up the rotation and torque direction is 
perpendicular to Bo and a positive torque induces a rotation 
of a negative 6. 



pie we investigate is a 2 mg bulk crystal, consisting 
in an ensemble of independent, iso-oriented and identi- 
cal rings. In fact, the synthesis procedures developed 
in supramolecular chemistry allow the engineering of 
regular arrays, as required (together with an interclus- 
ter coupling) by the long-term goal of scaling up the 
quantum hardware. The experimental details and the 
method are described in Ref. 0, while the synthesis 
and a systematic study of the thermodynamic properties 
of [{Me2NH2}{Cr7MF8(02CCMe3)i6}] (in short CryM, 
with M=Fe, Co, Ni, Mn, Cd) will be reported elsewhere. 
In Fig. n we show C{T;Bo), normalized to the gas con- 
stant R, for five different values of Bq. The typical low- 
T Schottky anomaly appears upon application of the 
magnetic field and is shifted towards higher tempera- 
tures as Bq increases. The overall specific heat C(T; Sq) 
arises from two distinct contributions: one is a lattice- 
related Debye term, Ciatt/R = 2MrT^/{QD + XT^f, 
where r = 298 is the number of atoms per molecule, 
Gd = (158 ± 10) K and A = (0.42±0.02) K"!; the other. 
Cm, depends on the spin degrees of freedom and is re- 
sponsible for the Schottky anomaly. The dependence of 
Cm on the energy eigenvalues of H, reads: 



Cm — R0^ 



S,:e" 



(2) 



with /3 = I/UbT. The best fit of the specific-heat tem- 



perature dependence at different values of Bq (solid lines 
in Fig. nj provides the following set of values for the mi- 
croscopic physical parameters: Jcr/ks — (17 i 0-5) K, 
Jm/Jcr = 0.9±0.1, \d,\/kB = (0.3±0.15) K (for all the 
spin sites) and^Ni = 2.2±0.1; gcr = 1.98±0. 02 was deter- 
mined from previous measurements on Crg. The sign of 
the anisotropics di needs to be independently determined, 
as we do by means of torque magnetometry measure- 
ments. In the inset of Fig. ^ we report an angular scan 
of the torque signal, measured at T = 2 K and i?o = 2 T. 
The torque exerted, e.g., for slightly off-perpendicular 
magnetic- field orientations {6 ~ 90°) tends to align the 
ring plane to the direction of Bg. For symmetry reasons, 
the doublet ground state gives no contribution to the 
torque, and those of the first excited multiplet (5 = 3/2, 
see below) are the only excited states to be meaning- 
fully populated at T = 2 K: the above behaviour there- 
fore shows that the S=3/2 manifold has an easy-plane 
anisotropy, which implies < 0. Results independently 
obtained by high-field torque measurements and inelastic 
neutron scattering (INS) confirm the above estimates of 
the microscopic parameters Lower-symmetry terms 
probably need to be included in the spin Hamiltonian in 
order to refine the interpretation of the INS data, but are 
neglected in the following since, they have little relevance 
to this discussion. Finally, the intracluster dipole-dipole 
interaction has been evaluated by means of the point- 
dipole approximation. 

We are now able to include the microscopic parameters 
in Ti. and accordingly draw the pattern of the low-lying 
energy levels as a function of Bq = BqZ (Fig. ^ . At 
zero field the ground state is a degenerate doublet (eo.i), 
with a largely dominating (> 99%) total-spin component 
S = 1/2; the first excited states (£2-5), instead, belong 
to a typical rotational band, with S ~ 3/2. Noticeably 
enough, these low values of the ^-mixing allow us to con- 
sider the total spin as a good quantum number for the 
lowest eigenstates of the CrrNi molecule. As discussed 
in more detail in the following, two other physical quan- 
tities also play a crucial role in the perspective of a QIP 
implementation. The first one is the energy difference 
between the ground-state doublet and the higher-lying 
levels, i.e. A = £2 — Ci: in fact, A determines to which 
extent the ring behaves as an effective two-level system, 
i.e. a meaningful population of any state but |0) and 
|1) can be avoided throughout the molecule manipula- 
tion. The second one is the splitting between the two 
S" = 1/2 states, 6 = ei — eq: S fixes the temperature 
the systems has to be cooled at in order for it to be 
initialized to its ground state. Likewise the same ex- 
istence of a ground-state doublet and the suppression 
of the ^-mixing, the large energy separation from the 
higher states in the present molecule, A(0) ~ 13 K, 
is a non-trivial result of the system engineering. Be- 
sides, the magnetic field allows a further tuning of the 
molecule's level structure. In particular, it increases S, 
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FIG. 2: Energy levels of the CryNi molecule as a function 
of a static magnetic field applied along the 2;-axis. At zero 
field the ground-state doublet is energetically separated from 
the higher states (A(0) ~ f3 K), thus representing a suitable 
choice for the qubit encoding. 

whereas it decreases the energy difference A{Bo) between 
|1) ~ IS* = 1/2, M = 1/2) and |2) ~ |3/2,-3/2). The 
achievement of the best trade-off between the conflicting 
requirements of maximizing d and A therefore determines 
the optimal value of the field, which we identify with 
Bo = 2 T (see the discussion below). As a consequence, 
the achievable temperature required for the system ini- 
tialization to the |0) state is T << S/ks — 2.4 K, whereas 
A{Bo) ~ A{0)-2gnBBo ^ 9.4 K. On the grounds of the 
above results, the CryNi molecule can be considered as 
an effective S — 1/2 spin cluster, and the information 
states |0) and |1) safely identified with its ground-state 
doublet. Such conclusion is unaffected by the uncertainty 
on A(0) (resulting from that on the microscopic param- 
eters) which is roughly equal to that on Jcr uM- 

The time simulation of the quantum gates pro- 
vides an important feedback for the optimization of 
the physical parameters. In fact, the general uni- 
tary transformation applied to the computational space 
is decomposed into a sequence of elementary gates, 
such as the SU{2) rotations of the single qubit and 
the two-qubit CNOT 11]. We start by considering 
the former, which can be obtained as a combination 
of 3 rotations about any two orthogonal axes, e.g. 
t/(a,/3,7) — exp(—iQ!cr2) exp(—i/3(Ti) exp(— 170-2), being 
ci,2,3 the Pauli matrixes. Transitions between the |0) 
and |1) states, i.e. rotations about the x and y axes, can 
be induced by means of resonant, in-plane electromag- 
netic pulses Bi{t) cosojt, where Bi{t) <C Bq represents 
the slowly-varying envelope. In the case of an effective 
two-level system, however, the transverse magnetic field 
also couples the |0) and |1) states to the higher-lying ones, 
thus inducing a population loss (leakage) during gating, 
quantified by L = 1 - [\{0\'ip{t))\^ + |(l|?/'(t))P ]. More 
specifically, the occurrence of such unwanted transitions 
is due to the S'-mixing and to possible intraclustcr in- 
homogeneities in the magnetic fields (or, equivalently, in 
the ion g factors): both result in non- vanishing matrix 

elements (0, l|7Yop|i > 2), with Hop = Ms J2i=i 9^ ^1 ' ^i- 
Together with a molecular engineering aimed at the sup- 
pression of the S'-mixing, the minimization of L can be 
achieved by the use of "soft" enough pulses, i.e. by keep- 
ing the pulse spectral dispersion Auj < {A~ d)/h. Being 
Alo ^/Tg, such inequality gives a lower bound to the 
the ratio Tg/rd that can in principle be achieved (t^.^ are 
the gating and decoherence time, respectively). 

Here, by means of a numerical integration of the 
Schrodinger equation, we calculate the time dependence 



of |co,ip = 1(0, l|V'(t))|^ and of L corresponding to a tt 
rotation about the x axis, for \tp{0)) = |0) and i3o = 2 T 
(see Fig. |3K); the pulsed field is assumed to have a 
gaussian profile Bi{t) = Aexp[—{t — to)^ /{2a^)], with 
A = 0.1 T and a — 150 ps. To a very large extent the 
leakage involves the first excited multiplet, and its value 
remains lower than 10~^ throughout the time evolution: 
the lower limit of Tg imposed by the presence of upper 
levels is thus by far larger than the one arising, e.g., from 
the present technological limitations in the pulse gener- 
ation. Further simulations, performed with more rea- 
sonable values of the mag netic field [A = 0.01 T) [l^, 
give proportional increases of the pulse duration, witht a 
further suppression of L. We have also investigated the 
gate robustness with respect to possible spatial inhomo- 
geneities of Bi, generally leading to a larger values of L: 
even in the worst limiting case of a Bi which is nonzero 
only at the Ni-ion site, L remains below the threshold of 
10-4 = 0.15 T, Tg ~ 2cr - 300ps). A larger value of 
the static field Bq would polarize nuclear spins, thus re- 
ducing the decoherence due to hyperfine field fluctuations 
(see below). For Bq > 5.6T, ei > £2 and |1) can relax into 
1 2) by emitting a phonon (the decay to |0) is practically 
forbidden, owing to the approximate Kramers-doublet 
nature of the computational basis). These processes, 
however, are expected to be completely negligible with 
respect to other sources of decoherence, due to the small 
value of (ei — €2) /ks^D- An increase of L might arise 
from the breakdown of the condition that A > Aw. In 
order to estimate this effect, we have considered a static 
field Bo = 5.61 T, which implies A ~ (Fig. Eb): L is 
seen to reach ~ 3%, but even small deviations of Bq from 
this critical values is enough to suppress L. 

The experimental demonstration of the intercluster 
coupling required for the implementation of the two-qubit 
gates is beyond the scope of the present work. In the 
following we discuss the present scenario and possible 
strategies to be pursued in order to achieve the required 
conditional dynamics. In particular, links between CryNi 
rings formed by delocalized aromatic molecules have al- 
ready been synthesized 13], and intercluster couplings 
have been demonstrated in similar systems 14]. The 
simplest case occurs when the coupling between the rings 
results from the interaction of the m-th spin in molecule 
/ with the n-th of //. If we apply the first-order per- 
turbation theory and restrict ourselves to the product 
space {|0), Jl)}/ ® {JO), Jl)}//, a Heisenberg interaction 
Hiji — J*sl^ ■ s^^ can be rephrased as follows in terms 
of the cluster total spins: 

HlJI = Jeff^I ■ S//, (3) 

where 

Jeff = 2/3 r (l/2lls„lll/2) (l/2lls„lll/2) (4) 

(analogous equations apply to the case of an Ising in- 
teraction). In the present molecules the single-cluster 
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FIG. 3; Simulated time evolution of the molecular magnet, 
initially prepared in |0), under the effect of a transverse mag- 
netic field with a gaussian temporal profile. The black, red 
and blue lines correspond to |coP, |cip and L (multiplied by a 
factor 10^ in (a)), respectively; the static field is (a) Bq — 2T 
and (b) Bo = 5.61 T, whereas ^ = 0.1 T in both cases. 

reduced matrix elements (l/2||sm,n||l/2) has a modulus 
of the order of unity and alternating positive-negative 
signs as a function of n and to, whereas it is always neg- 
ative for the Ni ion. If the interchain couplings are more 
than one, Jeff is given by the sum of the single contribu- 
tions. Therefore, Jeff has the same order of magnitude 
of the coupling J* between single spins, whereas its site 
dependence provides additional flexibility to the imple- 
mentation scheme. 

Within the present scenario the inter-cluster couplings, 
though engineerable during the growth process, are un- 
tuneable thereafter. This shortcoming is common to 
other implementations, and doesn't in principle pre- 
vent from performing the single- and two-qubit quan- 
tum gates. In fact, depending on the form and magni- 
tude of the interchain coupling, different approaches can 
be adopted. Weak and Ising-like interactions allow the 
use of the so-called "refocusing techniques", widely de- 
veloped within the NMR community Off-diagonal 
(e.g., Heisenberg) intercluster interactions could instead 
require either multi-ring encodings of the qubit, result- 
ing in a symmetry- induced cancellation of the logical cou- 
pling between the encoded qubits , or the use of inter- 
qubit molecules acting as tuneable barriers p^ . 

We finally discuss the decoherence of the cluster-spin 
degrees of freedom, which is expected to mainly arise 
from the hyperfine coupling with the nuclear spins. A 
first estimate can be obtained by considering the dipo- 
lar interaction of one Cr ion (s=3/2) with the neigh- 
boring F nucleus (natural abundance ~ 100% of the 
7=1/2 isotope). Being gp = -f-5.2577 and the dis- 
tance of each of the eight F atoms from the nearest Cr 
ion d = 0.191nm, the interaction energy corresponds to 
Ehyp/ks = 0.38 mK. For an octanuclear ring this would 
give Td ~ h/8Ehyp ~ 2.5 ns, that can be considered as 
a lower bond for r^. Similar values have been esti- 
mated for other molecular magnets 0, . Direct mea- 
surements of the electron-spin resonance linewidth on a 
CryNi crystal (which also includes the effects of possi- 
ble inhomogeneities, the so-called "dephasing") provides 
Td values one order of magnitude larger than the above 
one 19]. Substantial enhancements of will result 
from the suppression of the hyperfine field fluctuations 
(mK temperatures and few-T static fields) and from the 
substituting of the F ions with OH groups in the Cr7Ni 
compound, which we are already working at. Besides, 
the static field induces a large mismatch between the en- 



ergy gaps of the nuclear and electronic spins, thus ren- 
dering highly inefficient the relaxation processes of the 
latter. A second potential source of decoherence, namely 
the ring-ring dipolar coupling, is characterized by a lower 
energy scale, Ed^p/kB ^ {gtisSf/V - 0.1 mK {S = 1/2 
and V — 6.346nm'^), futher reduceable in diluted Cr7Ni 
molecular systems. 

In conclusion, the energy spectrum of the investigated 
CryNi molecule fully justifies its description in terms of 
an effective two-level system; besides, the symmetries 
of the ground-state doublet (S'-mixing below 1%) sup- 
press the coupling to the higher levels as induced by 
the transverse magnetic fields which are required for the 
quantum-gate implementation. In fact, our simulations 
of the single-qubit gates provide negligible values for the 
leakage [L < 10^*) even for gating times of the order of 
10^ ps, i.e. well below the tens of ns estimated for the 
spin decoherence times. While further work is needed for 
the engineering of the intercluster coupling, these results 
strongly support the suitability of the CryNi rings for the 
qubit encoding and manipulation. 
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